How long does it take to pull an ideal polymer into a small hole? 
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We present scaling estimates for characteristic times run and Thi of pulling ideal linear and ran- 
domly branched polymers of N monomers into a small hole by a force /. We show that the absorbtion 
process develops as sequential straightening of folds of the initial polymer configuration. By esti- 
mating the typical size of the fold involved into the motion, we arrive at the following predictions: 
Tnn(A^) ~ N 3 / 2 j f and r^ r (N) ~ N 5 ^ 4 /f, and we also confirm them by the molecular dynamics 
experiment. 
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There are models in physics - most famously exempli- 
fied by an Ising model - whose importance is not due to 
their practical applicability to anything in particular, but 
because they provide the much needed training ground 
for our intuition and for the development of our theoret- 
ical methods. In polymer physics, one of the very few 
such models is that of ideal polymer absorbtion into a 
point-like potential well. For a homopolymer in equilib- 
rium, this is a clean yet interesting example of a second 
order phase transition 0-Q For a heteropolymer, even 
in equilibrium, the problem is still not completely under- 
stood - see, for example, 0, El IE ; mathematically it 
is similar to the localization (pinning) transition in the 
solid-on-solid models with quenched impurities 0, ■ 

Surprisingly enough, the dynamics of a polymer chain 
pulled into a potential well is almost not discussed at all. 
On the first glance, one might think that this dynamics 
should be related to the dynamics of coil-globule collapse 
[ToL fTTl ll2L IT3I IbH ] , but the similarity between these two 
problems is limited. Indeed, globule formation after an 
abrupt solvent quench begins by a simultaneous chain 
condensation on many independent nucleation centers, 
while the absorbtion process can be viewed as pulling of 
a polymer with a constant force into a single hole. Hence, 
only a single attractive center causes the chain condensa- 
tion. In this sense, dynamics of absorbtion into a poten- 
tial well is more similar to a driven polymer translocation 
through a membrane channel 0, [n| EB [H3| • 

To describe in words the polymer pulled into a hole, let 
us imagine a long rope randomly dropped at the desk top. 
Pull now the rope by its end down from the desk edge. It 
is obvious, and can be established through an experiment 
accessible even to theorists, that the rope does not move 
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all at once. What happens is only a part of the rope 
moves at the beginning, such that a fold closest to the end 
straightens out. As soon as the first fold is straightened, 
it then transmits the force to the next fold, which starts 
moving, and so on. The rope sequentially straightens its 
folds, such that at any moment only part of the rope 
about the size of a current fold is involved in the motion, 
while the rest of the rope remains immobile. 

In the present letter we show how this mechanism ap- 
plies to an ideal homopolymer dynamics and how it sheds 
light on the physics behind scaling estimates 0] for 
the relaxation times for both linear, rn n , and randomly 
branched, Tb r , polymers. 

We address the problem on the simplest level of an 
ideal Rouse homopolymer of N monomers, each of size 
a, embedded into an immobile solvent with the viscos- 
ity rj. One end of the chain is put in the potential well 
(the "hole"), and the polymer is pulled down to the hole 
with constant force /, acting locally on just one single 
monomer, the one currently entering the hole. We will 
be interested in the relaxation time r of the complete 
pulling down of the entire polymer into the hole, we want 
to know how r depends on polymer length N and the 
pulling force /. Neglecting inertia, the balance of forces 
reads / = ff T (t), where friction force /f r should be linear 
in velocity with n(t) being the number of monomers 
'swallowed by the hole' by the time t. More subtly, fric- 
tion coefficient must be proportional to the number of 
monomers m(t) currently involved in the forced motion: 
/fr = T) a 2 m(t) dn Jp ■ Our main task will be to estimate 
the length m(t) involved in the motion for either linear 
or branched polymers. 

In order to clarify our approach, let us compute the 
characteristic time r of pulling down of an initially 
stretched linear polymer. In this case, all monomers 
move except those already absorbed in the hole, so 
m(t) = N — n(t), and so the force balance condition 
reads 



f = V a 2 [N-n(t)} 



dn(t) 
dt 
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Upon integration, this yields 
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In terms of dependence on N, this result is similar to the 
Rouse relaxation time of a linear polymer, tr ~ tqN 2 , 
where r ~ r/a 3 /ksT, and T is temperature. Our es- 
timate, (J2J, remains valid as long as r s traight < tr, be- 
cause in this case the initially straight polymer has no 
time to coil up on itself while being pulled up by the 
force /. Given the expression for r st raight, ©, the condi- 
tion r s traight < T R translates into / > fc^T/a, the latter 
making an obvious sense: pulling force should be strong 
enough to cause significant stretching of every bond, or, 
in other words, Pincus blobs must be as small as of the 
order of one monomer 20] . 

Now we turn to more realistic dynamics starting from 
a Gaussian coil configuration, and look at Tn n (7V). To 
begin with, let us consider a simple scaling estimate . 
First of all, we note that the only time scale relevant 
for the problem is the above mentioned Rouse time tr, 
which is the time needed for the coil to diffuse over the 
distance about its own size, R ~ aiV 1 / 2 . With tr be- 
ing the solitary relevant time scale, the absorbtion time 
should be written in the scaling form 
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where the factor 4> can depend on all other rele- 
vant parameters only in the dimensionless combination 
/i?/(fcsT). The second step follows from the fact that 
the speed of the process, N/m n , must be linear in the ap- 
plied force, which means that the function 4>(x) should 
be inversely proportional to its argument, <j){x) ~ 1/x. 
This then yields 
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Scaling argument gives the answer, but does not give 
much insight. To gain proper insight, let us show how 
this same answer results from the process of sequential 
fold straightening. The definition of "folds" and of their 
lengths Sf \(N) is illustrated in Fig. ^ Since the chain 
is being pulled all the time into an immobile point in 
the center, the definition of a fold arises from considering 
only the radial component of the random walk, |r(s)|, 
where r(s) stands for the (natural) parametric represen- 
tation of a path. In 3D, for example, |r(s)| represents 
Bessel random process jU HH- To obtain lengths of 



folds 
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Sfoi , one has to coarse grain the tra- 
jectory up to the one monomer scale a, thus making a 
Wiener sausage, and then find the local minima of |r(s)| 
separated by the biggest maxima. It is obvious that there 
are about y/N of such minima, separated by the intervals 
about y/N each. 

When we pull down the Gaussian chain by its end with 
a constant force /, only the current fold, of the length 
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1: (Color online) Random coil viewed as a sequence of 



N, is involved in the forced motion and expe- 
riences friction. The equation is still valid, but N 

should be replaced by s[*| and n should be integrated in 

yielding relaxation time for one fold 
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All ~ ViV folds relax sequentially, meaning that Ty m ~ 
\/~N Tfoi , and thus returning the scaling answer . 

Notice that the absorbtion time of a linear chain, Ti; n , 
Eq. Q, is much shorter than the typical Rouse time 
under the very weak condition / > kBT/(ay/~N), which 
means Pincus blob dictated by the force / should only be 
smaller than the entire coil. Under this condition, dur- 
ing the absorbtion process the configuration of the chain 
parts not involved into the straightening of a current fold, 
can be considered as quenched. 




FIG. 2: (Color online) Randomly branched chain and the 
corresponding hierarchy of scales. 

To understand the chain pulling dynamics even bet- 
ter, it is useful to consider another fundamental fractal 
model, that of randomly branched polymer, schemati- 
cally shown in FigEI T° nn d scaling approach for the ran- 
domly branched polymer, we have to estimate the analog 
of Rouse time as the time for a polymer to diffuse to its 
own size. Since all monomers experience friction indepen- 
dently of each other (there are neither solvent dynamics, 
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nor hydrodynamic interactions), it follows that the diffu- 
sion coefficient of a Rouse coil, linear or branched, is given 
by D r~j k B T '/ '(rjaN) . However, the radius of a branched 
polymer, unlike a linear one, is only R ~ aN 1 / 4 . There- 
fore, the Rouse relaxation time for the branched polymer 
reads tr ~ tqN 3 / 2 . The rest of the argument repeats 
exactly what we have done for the linear chain, yielding 
the result 
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Let us show now how one can re-derive (JHJ generaliz- 
ing the concept of Bessel process to randomly branched 
chains. It is convenient to introduce a hierarchy of scales 
in the iV— link ideal randomly branched chain. Let us 
denote by N\ a number of monomers in a typical " prim- 
itive" (or "bare") path of a largest scale. The size R of 
an ideal randomly branched polymer is R ~ aN 1 / 4 . The 
bare path of the maximal scale is a Gaussian random walk 
with R ~ aN^ 2 . Hence, N\ ~ N 1 / 2 . The average num- 
ber of monomers in all side branches N%, AT 3 , JV4, etc. 
attached to one bond of the bare path of the maximal 
scale (as shown in Fig|2Ji), is of order of N/Ni ~ N 1 / 2 . 
All these side branches form a randomly branched con- 
figuration of size R 2 ~ a(iV/A^i) 1 / 4 ~ aN 1 / 8 . Now, 
the number of monomers in the bare path of the second 
scale, N 2 , can be obtained repeating the above scaling 

arguments: R 2 ~ aN^ 2 , giving N 2 ~ N 1 / 4 , and so on: 
N 3 ~ iV 1 / 8 , . . . , N k ~ N 2 ~ k see FigH 

On all scales the bare paths Ni, N 2 , . . . Nk, etc. form 
the Gaussian sub-coils, each of which can be divided in 
folds Sfoj , Sfo? , . . . Sfo) , . . . in the same manner as it is 



(k) 

done above for a linear chain. For s fo ( we have the esti- 
mate 
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If the randomly branched polymer is pulled by its end 
with a constant force, /, then all Gaussian sub-coils on 
all scales straighten their folds simultaneously. The to- 
tal typical number of monomers s a u in all folds simul- 
taneously involved into the motion on all scales can be 
estimated as 
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^1/4+1/8+1/16+... = aN l/2 (g) 



We have seen on the example of the linear chain, that 
the characteristic time to straighten the fold of a typical 
length s is of order of r ~ s 2 - see ©. So, as it follows 
from (JSJ, the characteristic time of straightening out all 
folds of length s a n is of the order of 
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In a randomly branched polymer there are N\a/sQ 
independent parts, which relax sequentially. Hence, 



the total relaxation time can be estimated as Tb r 
s foi) r foi' which returns the result (JBJ. 



We have tested our results by a molecular dynamics 
experiment. We used standard Rouse model with har- 
monic bonds between neighboring monomers along the 
chain and with the averaged bond length, a = 1. The at- 
tracting force / acts at any moment on one monomer only 
and always points to the origin. Once a given monomer 
k gets inside the absorbing hole (|rfc| < R a = 1, where 
is the position vector of the kth monomer) this monomer 
is considered absorbed, it does not move any more and 
does not exert any force on monomer k + 1, its neigh- 
bor along the chain - but instead neighbor becomes the 
subject of the steady pulling force /. The same mecha- 
nism was simulated for the randomly branched polymer, 
except different branches can be pulled in parallel. The 
solution of dynamical equations is realized in the frame- 
works of the numerical velocity Verlet scheme [24j . The 
initial state of the linear chain was prepared by placing 
the first monomer at |ri| = R a . The initial configuration 
of the randomly branched polymer was prepared as de- 
scribed in the work |25| , and then polymer was shifted as 
a whole to place one randomly chosen monomer on the 
absorbing surface |r| = R a = 1. By performing 100 runs, 
we have computed averaged pulling times of linear chains 
from initially elongated and initially Gaussian conforma- 
tions, as well as for a randomly branched polymer. The 
results are shown in the Fig. [3]and demonstrate very sat- 
isfactory agreement with the theoretical predictions Q 
and iJOJ- 

In order to check more directly our proposed mecha- 
nism of absorption, we have looked at how Gaussian poly- 
mer chain gets involved into the biased motion driven by 
the pulling force /. The result is shown in the inset of 
figure |3 In this inset, k(t), marks the crossover between 
current fold and the rest of the chain. The current fold 
is operationally defined to include all monomers, with 
numbers j < k(t), whose motion is currently strongly af- 
fected by the pulling force /. Other monomers with the 
numbers j > k(t) are not yet affected. The simulation 
indicates that thus defined k(t) follows quite accurately 
the power law k(t) ~ i 2 / 3 with the exponent indepen- 
dent of the applied force. Inverting k(t), we get t ~ 
which suggests that the time of complete absorbtion of a 
/c-length subchain, which is given by our theoretical pre- 
diction and the time of involving fc-length subchain 
in the force-biased motion scale in the same way. 

To conclude, we have considered scaling estimates of 
the relaxation time associated with pulling linear or ran- 
domly branched polymer chain into a hole by a con- 
stant force /. We found that these estimates, rn n ~ 
(j]o? j f)N 3 / 2 and Tb r ~ (r) a 2 / f)N 5 / 4 , are consistent 
with the molecular dynamics simulation and they can be 
understood in terms of the scenario of sequential straight- 
ening of the polymer parts which we called folds and 
which were defined in terms of the radial component of 
the random walk representing the polymer chain. 

We have to emphasize that all consideration in this pa- 



4 




Elongated chain (el) 
Rand, coil (coil) 
Rand. br. polymer (br) 



10 100 
chain length, N 



FIG. 3: (Color online). Log- log plot of relaxation time against 
N. The units: a = 1, k B T = 1, and r\ = 10. Bullets (•) - ini- 
tially straightened linear chain, squares (■) - initially coiled 
linear polymer, triangles (a) - randomly branched polymers. 
The slopes are consistent with theoretical predictions. Inset: 
time dependence of the k(t) in the log-log scale. Here k(t) is 
the monomer number such that at the given time moment t 
monomers before k are already strongly involved in the biased 
motion towards the absorbing hole, while monomers after k 
are still diffusing not influenced by the force. 



per was restricted to polymer models without excluded 



volume in 3D. There is no doubt that self-avoidance, as 
well as hydrodynamic interactions will significantly af- 
fect the specific scaling relations obtained here. The ac- 
count of excluded volume effects is, to some extent, sim- 
pler than that of hydrodynamic interactions. On the ba- 
sis of scaling approach we can predict the pulling time, 
t(N) ~ N 1+v / f , for the polymer characterized by the 
exponent v in the relation R ~ N u . Our consideration is 
also restricted in the sense that we did not consider the 
role of the polymer part already 'swallowed' by the po- 
tential well. Although this might affect the result under 
some conditions, our aim was to elucidate the mechanism 
of sequential straightening of folds, for that purpose the 
absorbed tail is irrelevant. Similarly, in reality polymer 
is usually pulled into a hole in a wall, such as membrane, 
and we have neglected the (presumably logarithmic) fac- 
tors associated with the excluded half-space. 

However, as we said in the beginning, the purpose of 
the model is to facilitate methods and intuition. In this 
sense, we think that the consideration of ideal polymer in 
this letter was fruitful, because the mechanism of sequen- 
tial straightening of folds obviously applies to a number 
of real physical situations, such as, e.g., DNA transloca- 
tion through the membrane pore driven by the difference 
in chemical potentials. 
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